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FOR  INCOMPRESSIBLE  FLOWS 
By  C,  C.  Lin  and  S.  F.  Shen 


SUMMARY 


The  classical  theories  for  turbulent  shear  flow  are  the  momentum 
transfer  theory  of  Prandtl,  the  vorticity  transfer  theory  of  Taylor,  and 
®  the  similarity  theory  of  Von  Karman.  The  two  transfer  theories  both 

involve  a  mixture  length,  which  must  be  given  by  an  additional  assumption. 
^  On  the  other  hand,  the  similarity  theory  is  a  more  determinate  scheme, 

because  it  makes  a  more  definite  hypothesis  about  the  nature  of  the  turbu¬ 
lent  fluctuations.  Goldstein,  however,  introduced  an  alternative  form 
of  the  similarity  theory.  A  great  amount  of  work  has  been  done  to 
evaluate  the  relative  merits  of  these  three  theories. 

Further  investigation  into  the  nature  of  turbulent  motion  is,  however, 
done  largely  in  connection  with  the  simpler  case  of  isotropic  turbulence. 

In  this  field,  much  recent  progress  has  been  made,  particularly  following 
the  concept  of  Kolmogoroff.  The  concept  of  similarity  also  plays  a 
dominant  role.  Since  Kolmogoroff  *  s  theory  is  also  applicable  to  shear 
flow,  it  is  natural  that  one  should  reexamine  the  similarity  theory  of 
Von  Karman  by  using  modern  concepts.  This  is  the  main  purpose  of  the 
present  investigation.  It  is  found  that  the  original  form  of  the  theory 
is  supported  by  modern  concepts. 


INTRODUCTION 


The  concept  of  similarity  was  first  introduced  by  Von  Karman  in  1930 
^  (reference  l) .  Even  at  the  very  beginning,  he  realized  that  it  was  not 

possible  to  have  complete  similarity,  including  both  the  components  of 
%  fluctuation  essentially  free  from  viscous  forces  and  those  largely  influ¬ 

enced  by  viscous  forces.  It  was  Taylor  (reference  2),  however,  who 
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first  made  a  more  penetrating  analysis  into  this  question  in  connection 
with  isotropic  turbulence.  He  introduced  a  macroscale  of  turbulence  I 
and  a  microscale  X.  He  also  fully  discussed,  both  theoretically  and 
experimentally,  how  two  turbulent  fields  may  be  similar  in  all  the  large 
eddies,  which  contain  practically  all  the  energy,  and  yet  differ  com¬ 
pletely  in  their  rates  of  dissipation,  which  are  governed  chiefly  by 
the  small  eddies  having  a  negligible  contribution  to  the  total  energy. 

In  particular,  he  gave  the  relation 


(1) 


for  the  energy  dissipation  e  in  terms  of  the  scales  of  turbulence,  its 
intensity  u*,  and  the  kinematic-viscosity  coefficient  v.  This  rela¬ 
tion  shows  clearly  that  full  similarity  cannot  be  possible  in  general. 

A  further  advance  in  this  direction  was  made  by  Kolmogoroff  (refer¬ 
ence  3)*  According  to  his  concept^  for  large  Reynolds  numbers  of  turbu¬ 
lence,  as  defined  by  R^  =  u'X/v^  the  small  eddies  are  independent  of 

the  behavior  of  the  large  eddies,  except  to  the  extent  that  they  supply 
the  energy  to  be  dissipated.  There  are  then  only  the  two  parameters  6 
and  V  for  the  viscous  range.  From  dimensional  arguments,  Kolmogoroff 
introduced  the  characteristic  velocity  v  and  the  characteristic  scale 
r|  defined,  respectively,  by 


In  terms  of  the  spectrum  of  turbulence,  the  high  (spatial)  frequency 
components  are  dependent  only  on  v  and  T).  By  assuming  that  the  lower 
end  of  this  range  is  Independent  of  v  explicitly,  one  arrives  at  the 
spectrum 


(3) 


F  e2/3«-5/3 


where  k  is  the  wave  niomber.  This  relation  was  first  given  by  Obukhoff 
(reference  and  was  found  independently  later  by  Onsager  (reference  5)^ 
Heisenberg  (reference  6),  and  Von  Weizsacker  (reference  7)* 
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The  idea  that  the  low-frequency  range  takes  care  of  almost  all  of 
the  energy  and  that  the  high-frequency  range  takes  care  of  almost  all  of 
the  dissipation  has  been  fully  demonstrated  by  Von  Karmen  and  Lin  (refer¬ 
ences  8  and  9)*  It  is  found  that,  in  either  case,  the  unimportant  part 

-m  , 

to  the  important  part  is  of  the  order  of  (where  m  is  some  posi¬ 

tive  number)  and  is  consequently  negligible  for  large  Reynolds  numbers 
of  turbulence.  The  analysis  was  made  by  evaltiating  the  integrals 


=  F(k)  dK 


0 


I2  =  /  k:%(k:)  dK 

0 


J 


for  energy  and  dissipation,  assuming  the  characteristic  quantities  V,  2 
for  the  low-frequency  range  and  the  characteristic  quantities  v,  r\ 
for  the  high-frequency  range  and  using  Taylor’s  relation  (equation,  (l)) 
and  the  Obukoff  spectriim  formula  (equation  (3)). 

This  kind  of  analysis  is  now  used  to  remove  a  difficulty  raised  by 
Goldstein  in  his  analysis  of  Von  Karman’s  similarity  theory  (reference  lO) 
Goldstein  showed  that  there  are  at  least  two  ways  of  applying  the  simi¬ 
larity  theory,  one  analogous  to  the  momentum  transfer  theory  (the 
T-theory)  as  given  originally  by  Von  Karman,^  and  the  other  analogous  to 
the  vorticity  transfer  theory  (the  M-theory) .  By  an  analysis  of  the 
relative  importance  of  the  high-frequency  and  low-frequency  components, 
it  is  possible  to  show  that  the  r-theory  is  a  logical  consequence  of 
Von  Kerman’s  similarity  concept,  while  the  M-theory  does  not  follow 
directly.  This  kind  of  investigation  also  shows  that  the  usual  discus¬ 
sion  of  the  similarity  theory  needs  some  modification,  although  the 
final  conclusions  are  not  altered. 


It  is  the  purpose  of  the  present  paper  to  reexamine  critically 
Von  Karman’s  similarity  theory  for  incompressible  flows  by  using  modern 
concepts  in  order  to  provide  a  basis  for  extension  and  application  of  the 
theory  to  compressible  flows.  In  references  11  and  12  turbulent  boundary 
layer  over  a  flat  plate  in  compressible  flow  is  treated  in  the  same  spirit 
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This  investigation  was  carried  out  at  the  Massachusetts  Institute 
of  Technology  under  the  sponsorship  and  with  the  financial  assistance 
of  the  National  Advisory  Committee  for  Aeronautics. 

The  authors  wish  to  thank  Professor  S.  Goldstein  for  stimulating 
discussions  during  his  visit  here  in  September  1950. 

THE  CLASSICAL  SIMILARITY  THEORY  OF  VON  KARMM 


It  seems  convenient  to  begin  with  a  brief  sketch  of  the  original 
similarity  theory  of  Von  Karman,  together  with  Goldstein's  discussion 
of  its  difficulties  and  his  alternative  suggestion  (the  M-theory) . 
Remarks  will  then  be  made  from  the  present  point  of  view,  leading  to 
the  fuller  discussions  of  the  next  section. 

Von  Karmen  considered  a  steady  two-dimensional  mean  flow  and  made 
the  following  hypotheses  (see  reference  lO):  (l)  The  turbulence  mech¬ 
anism  is  essentially  independent  of  the  viscosity  of  the  fluid  (except 
in  the  viscous  layer  near  the  walls);  (2)  In  comparing  the  turbulence 
mechanisms  at  two  different  points,  consideration  of  the  fields  of 
turbulent  flow  may  be  restricted  to  the  immediate  neighborhoods  of  these 
points;  (3)  The  turbulence  flow  patterns  at  different  points  are  similar 
(relative  to  frames  of  reference  moving  with  the  mean  velocities  at  the 
points)  and  differ  only  in  scales  of  length  and  time  (or  velocity).  On 
the  basis  of  these  hypotheses  (the  validity  of  which  will  be  examined 
later),  the  following  development  of  the  theory  may  be  made.  (See 
appendix  A  for  definitions  of  important  symbols.) 


Consider  a  two-dimensional  parallel  mean  motion  with  velocity  U(y) 
in  the  direction  of  the  x-axis.  If  (u,v,w)  are  the  txirbulent  velocity 
components  and  p  is  the  pressxire,  the  Navier-Stokes  equations  are 


^  +  U  ^  +  l-(uu)  +  A(uv)  +  A(uw) 
St  Sx  Sx  Sy  Sz 


i  ^  +  V  Au 
P  Sx 


Sv 

St 


+ 


+  ^(vv) 

oy 


^(vw) 


1  ^  , 
P  hy 


V  Av 


^  (5) 


Sw 


+  u 


+  ^(wv) 

oy 


1  ^ 
P  Sz 


V  Aw 


fi 


IT 


NACA  TW  2541 


and  the  equation  of  continuity  Is 

^  ^  ^  =  0  (6) 

Sx  Sy  5z 

According  to  the  first  assumption  the  term  In  v  may  be  neglected.  One 
may  now  differentiate  the  set  of  eqiiations  (5)  and  combine  the  resultant 
equations  to  yield  the  equations  for  the  turbulence  vorticity  components 
usual  derivation  is  based  on  these  vorticity  equations. 

From  the  present  point  of  view,  the  neglect  of  viscosity  means  that 
the  ” low- frequency”  components  of  the  turbulent  fluctuations  are  being 
considered.  It  is  known,  however,  that  the  ’’high-frequency"  components 
of  the  vorticity  fluctuations  are  more  important.  It  is  consequently 
difficult  to  justify  the  neglect  of  viscosity  from  the  vorticity  equation, 
as  it  is  usually  done.  Even  the  neglect  of  viscosity  from  the  original 
Navier- Stokes  equations  (5)  requires  a  careful  examination.  These  points 
will  be  discussed  more  fully  in  the  next  section. 

*To  proceed  with  the  derivation,  take  the  origin  at  the  point  P 
under  consideration  and  axes  moving  with  the  mean  velocity  at  P,  so 
that  U  =  0  at  the  origin.  The  assumption  is  introduced  that  only  the 
immediate  neighborhood  of  P  may  be  considered:  For  dU/dy,  d^/dy^, 
their  values  are  taken  at  P,  while  for  U  the  first  term  only  of  the 
Taylor  expansion  is  taken,  namely,  ytdu/dy)^.  The  scales  of  length  and 

velocity  are  introduced  by  writing 


X  =  Zx 


t  =  Zt/A 


y  =  zy 


u  =  Au 


Z 


V  =  A  V 


so  that 


w  =  Zw 

I  =  A|/Z 
Ti  =  Arf/Z 

i  =  A?/Z 


where 


m 


and  so  forth.  Substituting  these  into  the  equations 
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for  vorticity  and  requiring  them  to  be  independent  of  the  position  of 
the  origin  and  of  the  values  of  dU/dy,  d^u/dy^  there ^  one  obtains 

A  =  Constant  I  — 

I  =  Constant 

Tlie  Reynolds  shearing  stress  at  different  points  of  the  fluid  is 

T  =  -puv  =  -pA^uv  and  is  therefore  proportional  to  pZ^(dU/dy)^;  and^ 
if  I  is  multiplied  by  a  suitable  constant^  there  may^  when  signs  are 
regarded,  be  written 


T  =  pl2 


dy  dy 


I 


y  du/d2u 
^  <iy/  dy2 


The  average  state  of  affairs  under  consideration  is  essentially  indepen¬ 
dent  of  X  and  z,  and  the  rate  M  at  which  x- momentum  is  communicated 
to  unit  volume  is 


M  =  d(-puv) /dy 


=  -ipA^/l)  d  (  uv  ^  dy 

and  is  therefore  proiwrtional  to  pZ(dU/dy)^  or  to  p2^(dU/dy)  (d%/dy^) . 

If  I  is  multiplied  by  a  suitable  constant,  there  may  therefore  be  written 


M  =  pZ2  ^  ^ 

dy  dy2 

I  =  Kg 


dU  /d^U 
^"/dy2 


A 
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Notice  also  that  M  =  if  the  average  state  of  affairs  is  indepen¬ 

dent  of  X.  Thus,  the  first  alternative,  the  T-theory  given  hy  Von  Karman, 
is  formally  related  to  the  momentum  transfer  theory,  while  the  second 
alternative,  the  M- theory  given  hy  Goldstein,  is  formally  related  to  the 
vorticity  transfer  theory,  although  the  basic  concept  of  the  similarity 
theory  is  essentially  different  from  a  transfer  theory. 

The  above  derivation  is  essentially  the  one  used  by  Goldstein,  the 
only  difference  being  that  he  started  with  the  vorticity  egiiations 
directly.  The  original  derivation  of  Von  Karmen  is  made  by  a  considera¬ 
tion  of  two-dimensional  fluctuations,  which  is  sufficient  to  provide  the 
essential  steps.  In  Goldstein’s  paper,  there  are  also  an  extension  of 
the  theory  to  the  flow  through  a  pipe  and  a  critical  discussion  of  the 
validity  of  the  theory.  His  investigations  show  that  in  the  axially 
symmetrical  case,  it  is  impossible  to  have  a  simple  similarity  as  that 
in  the  two-dimensional  case  if  terms  of  the  same  order  of  magnitude  are 
all  kept.  This  becomes  obvious  when  the  radius  of  curvature  of  the 
surface  of  constant  velocity  is  recognized  as  an  additional  parameter 
affecting  the  length  scale  for  a  local  similarity.  He  also  found  that 
the  velocity  distribution  from  the  t- theory  agrees  better  with  the 
experiments  for  flow  between  parallel  planes,  while  the  M- theory  appears 
to  be  more  satisfactory  for  pipe  flow. 

Another  generally  recognized  difficulty  with  the  similarity  theory  is 
the  following.  As  a  consequence  of  the  theory,  the  ratios 

u^:  v2:  uv:  uw:  vw  should  remain  constants.  In  a  channel,  this  is 

found  to  be  a  good  approximation  only  for  points  not  too  close  to  either 
the  wall  or  the  center  line. 

To  summarize,  the  following  unsatisfactory  points  in  the  similarity 
theory  have  been  discussed: 

(a)  There  is  nothing  in  the  original  theory  to  decide  between  the 
T-theory  and  the  M-theory  and  any  other  alternative  obtained^  say,  by 
an  application  of  the  theory  to  the  calculation  of  S^T/Sy^, 

(b)  It  cannot  be  extended  to  cases  other  than  the  original  case  of 
two-dimensional  parallel  flows. 

(c)  The  ratios  u^  :  v^  and  so  forth  are  not  constants  near  the  center 
of  the  channel  or  near  the  boundary, 

(d)  In  addition,  there  is  the  difficulty  in  applying  the  theory  to 
flows  with  a  point  of  inflection  in  the  velocity  profile. 

(e)  Again,  the  scale  I  turns  out  to  be  only  moderately  small  in 
most  cases. 
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It  must  te  noted  that  most  of  the  difficulties  are  limitations  which 
do  not  bear  on  the  basic  concept  of  the  theory.  The  final  form  of  the 
theory  is  applicable  only  when  the  scale  may  be  expected  to  be  determined 

by  dU/dy  and  d%/dy^.  This  is  obviously  not  the  case  in  wakes  or  in 
the  central  part  of  the  channel.  This  is  also  the  reason  why  the  axially 
symmetrical  case  cannot  he  treated.  The  presence  of  another  length  scale 
makes  similarity  unlikely.  Near  the  boundary,  similarity  might  also  break 
down  as  soon  as  the  scale  is  comparable  with  the  distance  from  the 
boundary. 

Thus,  the  essential  difficulties  are  (a)  and  (e)  discussed  above  if 
one  limits  oneself  to  flows  through  channels  and  in  the  boundary  layer. 

'It  is  to  be  recognized  that  the  whole  boundary  layer  or  the  channel 
should  be  regarded  as  an  "organic"  field  of  motion.  One  might  even 
describe  the  turbulent  motion  as  a  nonlinear  oscillation  superposed  over 
a  field  of  flow.  Thus,  any  attempt  to  "localize"  the  theory,  such  as  is 
done  in  the  similarity  theory  and  to  a  certain  extent  in  the  transfer 
theories^  is  at  best  a  rough  approximation.  One  is  faced  with  the  dilemma 
of  either  treating  each  individual  case  separately  with  proper  emphasis 
on  the  influence  of  the  boundary  conditions  or  being  satisfied  with  an 
approximate  theory  having  a  fairly  general  applicability.  The  latter 
course  has  often  been  taken. 

The  only  remaining  difficulty  (a)  is  to  be  settled  by  a  critical 
examination  of  the  concept  of  similarity  from  the  standpoint  of  recent 
developments  of  the  statistical  theory  of  turbulence.  It  will  be  seen 
that  the  basic  concept  of  Von  Karmfin’s  similarity  prefers  the  original 
form  of  the  theory  (t~ theory)  to  any  theory  based  on  higher-order  deriva¬ 
tives  (M- theory,  etc . ) . 


CRITICAL  DISCUSSION  OF  SIMILARITY  CONCEPT  FROM  STANDPOINT 
OF  MODERN  STATISTICAL  THEORIES 


The  concept  of  the  similarity  theory  in  shear  flow  will  now  be 
formulated  and  it  will  be  shown  how  it  is  related  to  the  concept  of 
similarity  developed  in  the  statistical  theory  of  isotropic  turbulence. 
The  concepts  developed  from  the  statistical  theory  will  then  be  applied 
to  the  present  problem  to  show  that  the  M-theory  does  not  follow  from 
Von  Karman*s  concept  of  nonviscous  similarity  in  shear  flow. 

General  concept  of  similarity.-  As  discussed  above,  the  similarity 
of  the  large-scale  eddies  (which  are  responsible  for  transfer)  is  at 
best  a  rough  approximation.  The  similarity  of  the  small  eddies  (which 


are  responsible  for  the  energy  dissipation)  is  however  a  much  closer 
approximation  at  very  large  Reynolds  ntimbers  of  tirrhulence.  According 
to  the  concept  of  Kolmogoroff,  this  part  of  the  turbulent  fluctviation 
will  depend  only  on  the  kinematic  viscosity  coefficient  V  and  the  rate 
of  energy  transferred  to  these  scales  from  larger  scales.  Furthermore, 
this  part  of  turbulent  motion  is  isotropic  and  has  a  loniversal  character 
independent,  say,  of  the  amoimt  of  shear  in  the  main  flow.  The  rate  of 
energy  transfer  is  approximately  the  same  as  the  total  rate  of  dissipa¬ 
tion  €.  TTm  order  of  accuracy  of  such  approximations  has  been  estimated 
by  Von  Karman  and  Lin  (references  8  and  9)*  In  other  words,  in  the  case 
of  shear  flow,  there  is  a  production  of  turbulent  ener^  from  main  motion 
by  the  usual  transfer  mechanism  (at  a  scale  of  Z  =  KU’/U",  say).  This 
large-scale  turbiilence  breaks  down  into  motions  at  smaller  scales  and 
eventiially  passes  into  the  viscous  range  to  be  dissipated  into  heat.  The 
rate  of  transfer  may  be  estimated,  say,  by  Taylor's  formula.  These  large- 
scale  motions,  being  produced  directly  from  the  mean  motion,  must  depend 
on  their  characteristics.  According  to  the  concept  of  Von  K^man,  the 
turb\ilent  fluctuations  in  this  range  have  a  universal  structure,  with  a 
length  scale  Z  =  KU'/u"  and  a  velocity  scale  ZU'.  Thus,  two  regimes 
of  similarity  are  visualized:  (a)  The  large-scale  similarity  of 
Von  Klirman  and  (b)  the  small-scale  similarity  of  Kolmogoroff.  The  first 
range  is  anisotropic  and  contributes  to  the  shear>  the  second  range  is 
isotropic  and  contributes  only  to  the  dissipation.  The  transition  range 
probably  has  an  energy  spectrum  62/3k"5/3  with  decreasing  amount  of 
shear. 

This  picture  of  similarity  is  analogous  to  the  one  visiialized  by 
Von  Kerman  and  Lin  (references  8  and  9)  for  the  intermediate  stage  of  decay 
of  isotropic  turbulence.  There,  the  large-scale  eddies  are  also  isotropic, 
having  a  scale  determined  by  the  Loltslansky  invariant.  It  must  be 
remarked  that  the  establishment  of  an  equilibrium  or  quasl-equllibrlum 
state  for  large-scale  motions  takes  a  long  time.  Thus,  it  may  be  sub¬ 
jected  to  doubt  whether  the  idealized  pictuire  thus  visualized  can  actually 
occur  for  decaying  turbulence.  On  the  other  hand,  in  the  case  of  shear 
flow,  a  stationary  system  being  considered,  the  condition  is  more  conducive 
to  the  establishment  of  Von  Karman' s  similarity. 

With  this  general  concept  of  similarity  in  mind,  the  classical  theory 
of  similarity  for  shear  flow  may  now  be  examined  with  the  help  of  results 
developed  for  isotropic  tvirbulence.  In  the  first  place,  one  may  examine 
the  orders  of  magnitude  of  length,  velocity,  and  vorticity  for  the  eddies 
of  large  scales  and  small  scales. 
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For  large-scale  motions, 


Length  scale 

1  - 

Velocity  scale 

Y 

Vorticity  scale 

v/z 

U'/U" 

lU' 


•scale  motions. 

Length  scale 

Velocity  scale 

V  -  - 

Vorticity  scale 

=  Nvimber  X  IR, 


-3/2 


y  (7) 


^  (8) 


where  is  the  Reynolds  number  of  turbulence  and  is  usually  very 

large . 


These  order -of -magnitude  relations  will  now  be  applied  to  the 
development  of  the  similarity  theoiy-  In  particular,  the  following 
two  points  will  be  considered: 

(a)  The  neglect  of  the  viscous  components  from  the  equations  for 

turbulent  fluctuations 

(b)  The  relative  plausibility  of  the  r-theory  and  the  M- theory 

In  making  these  considerations,  the  low-frequency  and  the  high-frequency 
components  will  be  considered  as  behaving  independently  in  a  linear 
equation.  For  nonlinear  terms,  low-frequency  components  can  be  obtained 
by  the  product  of  two  terms  both  of  the  high  frequency  or  both  of  the 
low  frequency,  while  high-frequency  components  can  be  obtained  by  the 
product  of  one  term  of  the  high  frequency  and  one  term  of  the  low  fre¬ 
quency.  Thus,  if  one  rites 

u  =  u^  +  u^  (9) 
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where  denotes  the  low-frequency  components,  and  the  high- 

frequency  components,  then  in 

the  low-frequency  components  come  from  u^^  and  u^^,  while  the  high- 
frequency  components  come  from  and  u^^. 

Since  u^  «  u^,  it  is  seen  that  u^  may  he  approximated  hy 

Thus,  in  the  left-hand  side  of  any  one  of  equations  (5)>  one  is  justified 
in  putting  a  subscript  I  to  every  fluctuating  quantity,  if  the  low- 
frequency  components  of  the  whole  left  side  are  desired.  The  right-hand 
side  is  linear,  and  hence  one  can  again  consider  the  low-frequency  com¬ 
ponents  separately.  For  these,  the  viscous  dissipation  is  known  to  be 

negligible  (as  may  be  verified  by  noting  that  V  Au  %  w/l  is  small 
compared  with  ^("au)  when  VZ/v  is  large).  Consequently,  Von  Karman*s 
original  assumption  of  neglecting  v  is  justified. 

The  equation  of  continuity  (6)  is  linear  and  can  therefore  be 
considered  separately  for  components  of  low  and  high  frequencies. 

It  is  remarked  that  the  above  discussion  can  be  carried  through 
only  when  equations  (5)  are  put  into  the  form  given.  If  the  left 
sides  of  equations  (5)  were  written  in  the  form 


Du 

Dt 


St 


^  +  U 


Su 

Sx 


(11) 


the  neglect  of  the  influence  of  the  high-frequency  components  on  the  low- 
frequency  components  would  have  been  dubious.  Consider,  for  example  the 
term  ^ 


V 


5y  hy 


Su- 


+  V 


h 


+  Vi 


ST 


(12) 
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m,  ,  ^ 

The  low-frequency  components  are  given  hy  vj  — i  as  well  as  hy  v,  _ _ 

^  ^y 

Now,  by  equations  (7)  and  (8), 


'h 


8u] 

8y 

8uj 

Sy 


YZ 

i 

T1 


I  ^ 


> 


(13) 


and  it  would  not  be  legitimate  to  replace  v  ^  by  v^^Su^/^y)  .  The 

mathematical  argiiment  can  he  carried  through  only  when  Du/Dt,  and  so 
forth  are  put  into  the  form  given  in  equations  (5),  with  the  help  of  the 
equation  of  continuity. 


T- theory  and  M- theory. -  The  above  arguments  also  suggest  that  in  the 
evaluation  of  r  =  -puv,  the  contribution  of  the  low-frequency  components 
predominates,  while  that  of  the  high-frequency  components  is  negligible. 
This  would  Justify  the  application  of  the  similarity  concept  to  the  cal¬ 
culation  of  T.  On  the  other  hand,  in  the  product  v^,  the  high- 
frequency  components  predominate.  This  suggests  that  would  be 

dependent  more  on  the  high-frequency  components  than  on  the  low-frequency 
components.  Thus,  the  similarity  concept  at  low  frequencies,  as  developed 
above,  cannot  be  used  for  obtaining  a  formula  for  The  M- theory  is 

therefore  not  a  direct  consequence  of  Von  Kerman’s  concept  of  similarity. 

These  arguments  should  be  somewhat  modified  by  the  fact  that  turbu¬ 
lence  tends  to  be  isotropic  at  high  frequencies.  This  fact  strengthens 
the  argument  for  t  but  weakens  the  argioment  against  v^.  It  is  pos¬ 
sible  that  the  contribution  to  v^  is  equally  important  from  high- 
frequency  and  low-frequency  components.  But  there  is  no  convincing 
argument  to  show  that  the  high-frequency  components  can  be  entirely 
neglected.  To  make  the  ideas  more  precise,  use  will  now  be  made  of  the 
correlation  tensor  and  the  spectral  tensor  in  the  discussion. 
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For  convenience,  first  consider  homogeneous  anisotropic  turbulence. 
Then,  the  correlation  tensor  =  Uj^(P)u^^(P')  is  a  function  of  the 

relative  positional  vector  The  spectral  tensor  stands 

in  Fourier  transform  relation  to 


If  a  spherical  coordinate  system  is  considered  in  the  K-space,  integration 
with  respect  to  the  angular  variables  leads  to 


F(k)cp^^(k)  dK 


(IT) 


where  is  a  "correlation  coefficient,"  being  equal  to  unity  when 

i  =  k.  Isotropy  at  high  frequencies  requires  that,  for  i  /  k, 
approach  zero  rapidly  with  k— »oo,  being  substantially  zero  in  the 
viscous  range. 


ThuS;  the  contrihution  of  the  low-frequency  components  to 

dominates  that  of  the  high-frequency  components  even  to  a  greater  extent 
when  i  ^  k  than  when  i  =  k.  In  the  latter  case^  it  is  well-known 
that  the  high-frequency  components  are  negligible  (cf*  appendix  B) ,  and 
this  is  therefore  still  more  so  in  the  case  i  ^  k  where  the  factor 
reduces  the  influence  of  the  high-frequency  components  still  further. 

It  is  not  possible  to  make  a  similar  discussion  for  by  using 

homogeneous  and  anisotropic  turbulence,  because  this  term  is  zero.  How¬ 
ever,  one  may  write  a  formula  analogous  to  equation  (l?)  in  the  form 


dK 


(18) 


where  is  a  dimensionless  "correlation  coefficient,"  approaching 

zero  rapidly  as  k — »«),  and  the  factor  k  is  suggested  by  formula  (15)* 

It  will  be  shown  in  appendix  B  that,  without  the  factor  the 

contribution  to  the  above  integral  would  practically  all  come  from  the 
high-frequency  components.  However,  since  ^±]^i  niust  approach  zero 
rapidly  as  k — ><»,  this  argument  ^  not  certain.  But,  at  any  rate, 
equation  (l8)  does  not  show  that  v^  will  be  essentially  determined  by 
the  low-frequency  components.  One  may  also  argue  as  follows.  The  fac- 

tor  probably  begins  to  become  insignificant  ohly  in  the  K  range 

For  various  Reynolds  numbers  of  turbulence,  corresponding  to  identical 

“5/3 

low-frequency  components,  the  extent  of  the  k  range  differs.  Since 

this  range  is  now  obviously  important  for  the  determination  of  inte¬ 
gral  (l8),  the  amount  of  vorticity  transfer  would  depend  on  the  local 
Reynolds  number  of  turbulence.  This  is  another  way  of  stating  that  there 
is  no  similarity  in  this  sense  between  the  various  points  of  the  flow 
field. 

The  above  discussion  shows  that  the  t- theory  should  be  used  in 
connection  with  Von  Karman’s  similarity  concept. 


DISCUSSION 


It  is  perhaps  in  order  now  to  make  a  general  survey  of  the  present 
theories  of  shear  flow.  In  contrast  with  the  similarity  theory,  there 
are  the  theories  of  transfer  of  momentum  and  vorticity.  In  comparing  the 
similarity  theory  with  transfer  theories,  it  should  be  noted  that,  although 
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there  is  formal  similarity  in  the  final  formulas^  there  is  a  definite 
difference  in  basic  concept.  The  similarity  theory  does  not  directly 
use  the  physical  picture  of  the  mechanism  of  transfer.  It^  however, 
asserts  more  definitely  on  the  nature  of  the  turbulent  fluctuations. 

On  the  one  hand,  this  has  the  advantage  of  leading  to  a  formula  for  the 
scale.  On  the  other  hand,  one  can  apply  the  concepts  evolved  from  modern 
statistical  theories  to  the  theory  of  shear  flow  by  following  the  concept 
of  similarity.  The  present  worh  seems  to  be  one  of  the  few  attempts  in 
this  direction,  and  it  tends  to  bear  out  the  classical  form  of  the  theory 
of  Von  Karman  (t- theory).  The  use  of  Von  Karman^s  concept  to  calculate 
average  quantities  involving  velocity  derivatives  (such  as  in  the 
M-theory)  cannot  be  justified  by  current  concepts  of  the  statistical 
theory  of  turbulence. 

The  conclusion  that  the  t~ theory  (analogous  to  the  momentum  transfer 
theory)  is  preferable  to  the  M- theory  (analogous  to  the  vorticity  transfer 
theory)  perhaps  requires  further  clarification.  A  great  deal  of  work 
has  been  done,  following  Taylor,  which  shows  that  the  vorticity  transfer 
theory  is  better  than  the  momentum  transfer  theory.  The  strongest  case 
is  perhaps  the  one  Involving  joint  velocity  and  temperature  distributions 
in  a  wake  and  in.  a  jet.  Whereas  the  momentum  transfer  theory  predicts 
the  same  distribution  for  velocity  and  temperature,  the  vorticity  transfer 
theory  can  account  for  the  difference  in  distribution  which  is  actually 
observed  experimentally.  In  fact,  if  0  is  the  fluctuation  of  tempera¬ 
ture,  the  transfer  of  heat  is 


.0v  =  2^ 


dy 


dy 


in  either  theory.  The  momentum  transfer  theory  gives  further 


T 

P 


-uv 


dU 


dy 


(19) 


(20) 


while  the  vorticity  transfer  theory  gives 


1  St  ^  ,2  dU  d% 


(21) 


If  I  is  assumed  constant,  this  leads  to 


dU 

dy 


(22) 
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The  difference  in  the  numerical  factor  1/2  is  the  essential  point  in 
question. 

€ 

In  applying  this  kind  of  argument  to  the  similarity  theory,  one 
must  recognize,  as  mentioned  above,  that  the  mechanism  of  transfer  is 
not  directly  used  in  the  similarity  theory.  In  the  transfer  theories, 
it  is  natural  to  use  the  same  transfer  coefficient  in  all 

the  formulas,  such  as  equations  (20)  to  (22).  On  the  basis  of  the 
similarity  concept. 


u  ~  Z 


V  ~ 


dy 


0 


i 


dy 


According  to  the  T-theory,  then 

dy 


dy 


v0  ~ 


dy  dy 


but  nothing  can  be  said  of  the  constants  of  proportionality.  Thus,  the 
above  method  of  testing  cannot  distinguish  between  the  r-theory  and  the 
M- theory. 

It  must  be  added  that  the  similarity  theory  with  the  for¬ 
mula  Z  ~  U’/u"  does  not  apply  to  the  case  of  wakes  and  Jets,  because 
the  velocity-distribution  curve  has  a  point  of  inflection.  The  general 
concept  of  similarity  of  low-frequency  components^  however^  might  still 
apply  even  though  the  scale  of  similarity  in  such  cases  is  probably  not 
detennined  by  the  local  velocity  distribution.  In  fact^  even  in  other 
cases  of  shear  flow^  the  validity  of  the  general  considerations  of 
similarity  does  not  depend  on  the  correctness  of  Von  Karmen's  formula. 


p 


w. 
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APPENDIX  A 

SYMBOLS 

I  similarity  scale  of  length 

p  pressure 

t  time 

u^v^w  fluctuating  velocity  components  in  x-,  y-^  and  z-direction^ 

respectively 

x^y^z  Cartesian  coordinates}  x-axis  in  direction  of  mean  flow 
L  macroscale  of  turbulence 

Reynolds  number  of  turbulence 
U  velocity  in  x- directi on 

V  velocity  scale 

^  rate  of  dissipation 

0  fluctuation  of  temperature 

K  wave  number 

X  Taylor’s  microscale  of  turbulence 

V  coefficient  of  kinematic  viscosity 
turbulence  vorticity  components 

p  density  of  fluid 

T  shearing  stress 

Subscripts: 

I  low-frequency  part  of  fluctuations 

h  high-frequency  part  of  fluctuations 

Barred  quantities  always  represent  mean  values;  primed  quantities 
represent  fluctuations.  - 
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APPENDIX  B 

MOMENTS  OF  THE  SPECTRUM  OF  TURBULENCp: 


As  indicated  in  the  section  "Critical  Dissussion  of  Similarity 
Concept  from  Standpoint  of  Modern  Statistical  Theories,"  the  relative 
contribution  of  the  low-frequency  and  high-frequency  components  to  the 
moments  of  the  spectral  function  plays  an  important  role  in  deciding 
the  physical  concept  of  the  mechanism  of  turbulence.  In  this  section, 
it  will  be  shown  that  in  the  case  of  large  Reynolds  numbers,  the  low- 
frequency  contribution  to  the  integral 


In  = 


F(k)k^ 


dK 


(Al) 


is  much  more  important  than  the  high-frequency  contribution  if  n  <  2/3, 
while  the  reverse  is  true  if  n  >  2/3  (equation  (A15))*  The  low  fre¬ 
quencies  and  the  high  frequencies  are  separated  at  a  fairly  arbitrary  r 

point  in  the  k  range.  In  the  borderline  case  of  n  =  2/3,  the 

relative  contribution  depends  on  where  the  separation  of  the  two  ranges  $ 

is  made  (equation  (A16) ) - 

The  special  cases  of  n  =  0  (giving  the  energy)  and  n  =  2  (giving 
the  dissipation)  are  well-known  (cf.  "Introduction")*  The  case  n  =  1  has 
been  used  in  the  section  "Critical  Discussion  of  Similarity  Concept  from 
Standpoint  of  Modern  Statistical  Theories"  to  bring  out  the  difficulties 
associated  with  the  M-theory.  The  case  n  =  -1  would  be  of  Importance 
in  determining  the  macroscale.  It  is  clear  that  the  macroscale  would  be 
proportional  to 


(A2) 


and  is  consequently  a  low-frequency  property.  This  is  in  contrast  with 
the  microscale  X  of  Taylor.  Since  X  is  proportional  to  I^/Ig,  it 
depends  on  the  low-frequency  components  as  well  as  on  the  high-frequency 
components. 

The  method  of  investigation  is  the  one  used  by  Von  Karman  and  Lin 
(references  8  and  9)*  some  frequency  (as  yet  unspecified)  in 
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the  range  vhere  F  ~  The  spectrum  helow  this  frequency  k*  ‘  is 

assumed  similar  with  a  length  scale  L  and  a  velocity  scale  V.  Above 
Itj  the  characteristic  quantities  are  v  and  q.  Thus^  the  spectrum 
may  be  described  as  follows.  For  low  frequencies, 


F(k)  =  ¥2Lf(KL)  (A3) 


where  the  function  f(x)  has  the  behavior 


-5/3 

f(X)  «  CX  '  (Al;) 


for  large  values  of  X,  for  example,  of  the  order  of  X*  =  k*L. 
For  high  frequencies. 


F(k)  =  v2T]g(KTl)  .  (A5) 


where  the  function  g(x)  has  the  behavior 


g(x) 


cx 


-5/3 


(a6) 


for  small  values  of  x,  for  example,  of  the  order  of  x*  =  K*q. 
In  evaluating  the  integral  (A1),  split  it  up  into  two  parts: 


P#  n  / 

In  “  /  F(k)k  dK  +  /  F(k)k^  dK 


(AT) 


The  low-frequency  part  is 


In,Z  =  /  7{k)k^  dK 


■r  ■■ 


f(x)  dX 


(A8) 
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where  X*  Is  larger  The  high-frequency  part  is 


F(k)k'^  dK 


=  —  I  x’^g(x)  dx 


where  x’*’  is  small.  To  get  the  order  of  magnitude  of  the  integrals, 
use  is  made  of  equations  (a4)  and  (a6)  .  Thus,  for  n  >2/3, 


x’^f(X)  dX 


x’^JfCx)  -  CX"5^  dX  +  CX*  - 


(AlO) 


The  first  integral  may  be  approximated  by  letting  the  upper  limit  go  to 
infinity,  and  then 


j  2  =  ax 


0(1) 

2 

n - 

3 


(All) 


In  the  case  n  <  2/3,  integral  (AlO)  is  convergent  as  X* — Thus 
write 


'  x’^f(X)  dX  =  I  x’^f(X)  dX 

0  Jo 


poo 

-  X’^  fl 

Jx* 


/  X  “5/3l  C  / 

f(X)  -  CX  dX  +  ^-(X*) 

J  n-- 
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which  leads  to  the  same  answer  (All) .  Similarly, 


x’^g(x)  dx  = 


2 

^  +  0(1)  for  n  ^2/3 

n  -  2 


The  ratio  of  the  two  integrals  2  '^n,h  therefore 


Jn,l/"n,h  -  op*' 


Jn.l/Jn.h  =  0  (i*) 


On  the  other  hand,  from  equations  (A3)  to  (a6) , 


2  /  x-5/3 

y‘^L  =  c(idL) 


v^c(kt)) 


n>| 


n<f 


SO  that 


Thus, 


■  (r 


,  r  -(n-lT 


n>| 


2 

n<- 


These  are  the  statements  made  at  the  beginning  of  this  section. 
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In  the  case  n  =  2/3,  the  above  type  of  arguments  give 


and  the  ratio  is 


In,l»  -573  1°SX» 


T1 


•n,l/^n,h 


log  r 

iiog  x*r 


(A16) 


which  depends  on  the  choice  of  the  frequency  k** 
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